Andreev reflections on Y\_ x Ca x Ba2Cu^0T-s evidence for an unusual proximity effect 
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We have measured Andreev reflections between an Au tip and Yi- x Ca x Ba2Cus,OT^s thin films 
in the in-plane orientation. The conductance spectra are best fitted with a pair potential having the 
d x 2_ y 2+is symmetry. We find that the amplitude of the is component is enhanced as the contact 
transparency is increased. This is an indication for an unusual proximity effect that modifies the 
pair potential in the superconductor near the surface with the normal metal. 

PACS numbers: 74.45.+C, 74.72.Bk 



The superconducting state is characterized by a pair 
potential (PP) A (r) = V (r) F(r) where V(r) is the po- 
tential representing the electron-electron interactions and 
F (r) = (r) ^| (r)) is the probability amplitude for 
finding a Cooper pair. In the case of a contact between a 
normal metal (N) and a superconducting metal (S), far 
from the interface we expect F(r) to be zero in N and to 
reach a constant value in S. However near the interface 
F(r) has a non zero value in the normal metal due to 
electron pairs leaking from S into N. This phenomenon 
is known as the proximity effect and has been studied 
both theoretically and experimentally for the case of su- 
perconductors with a pair potential having an s-wave 
symmetry [lj. Microscopically the mechanism involved in 
the creation of this non zero pair amplitude F(r) in the 
normal metal, is a special reflection process known as An- 
dreev reflection [2j, which occurs when there is an abrupt 
change in the PP. Such a change occurs at a clean N/S 
contact, where F(r) is continuous, due to the difference 
in V(r) in the two materials. An electron approaching 
the superconductor(SC) from the normal metal with en- 
ergy smaller than A cannot enter as a quasiparticle into 
the superconducting condensate. Instead the electron is 
reflected as a hole and a Cooper pair is added to the con- 
densate. Therefore by measuring Andreev reflections one 
can investigate the properties of the superconducting PP. 

For a perfectly transparent, small N/S contact, of size 
a << /, where I is the mean free path, the Andreev re- 
flection process is manifested by a low bias (eV < A) 
conductance which is twice as large in comparison to the 
high bias one. Blonder et a/, jjjj] (BTK) have calculated the 
conductance of an N/S contact with an additional barrier 
represented by a delta function potential U (x) — H5 (x) . 
They have defined a dimensionless parameter represent- 
ing the barrier strength, Z = ^^ r , where vp is the Fermi 
velocity. A clean N/S contact is described by Z=0 and 
a high barrier tunnelling contact by Z >> 1. For finite 
Z values, the conductance is depressed at low bias and 
enhanced at eV = A. In the case of high temperature su- 
perconductors (HTS) various experiments have indicated 
that the PP has a d x 2_ y 2 (d) symmetry 0, described by 
A (k) = Aq cos (28), where 9 is the polar angle measured 



from the crystallographic a-axis. This PP is very dif- 
ferent from the isotropic s-wave PP as it changes sign 
and has nodes for 6 = ^ + \ n, where n is an integer. 
Thus it is expected that the properties of N/HTS con- 
tacts should differ from those of N/S contacts. Tanaka 
and KashiwayajS] have extended the BTK calculation to 
the case of an anisotropic PP, A(8) and have shown that 
the conductance curves differ from those calculated by 
BTK. They have defined Z = = 2Z BTK hence- 

forth we shall use this definition. For a d PP a zero 
bias conductance peak (ZBCP) appears for in plane con- 
tacts reflecting the existence of Andreev surface bound 
statesQ. It is most pronounced for (110) contacts. For a 
(100) contact, the low bias dip which evolves in the case 
of a s-wave PP for Z > 0, appears only for Z > 0.4 and 
the amplitude of the normalized conductance maxima is 
lowered considerably, reaching a value of around 1.5, as 
shown in Fig. la, compared to a value of around 2.1 in 
the s-wave case for this Z value. Experimentally, data 
on Andreev reflection at low Z (Z<l)contacts with HTS 
is limited. Results on YBCO were reported by Hess et 
aL0 and Yagil et al. fl(i| . Wei et al. 11] have reported 
a measurement on YBCO with a fit to a d PP. Mea- 
surements on La2- £C Sr :! ;Cu04 were reported by Achsaf 
et aL|l2| with a fit to an anisotropic s-wave PP and by 
Gonnelli et aZ.^j| with a fit to a d+id xy PP. 

Here we report on Andreev reflection spectroscopy 
measurements using a point contact between a normal 
metal (Au) and Yi- x Ca x Ba2Cu^O'j-s thin films. By fit- 
ting our data to the theory of Tanaka and Kashiwya0, 
we are able to find the symmetry of the PP at different 
barrier transparencies. Our results are best fitted to the 
model using a PP having the d+is symmetry. The am- 
plitude of the is component is large at low Z (Z < 0.5), 
where it reaches 80% of the d amplitude, and small at 
large Z. This dependence on the barrier transparency im- 
plies that the large value of the is component seen at 
low Z values is the result of an unusual proximity effect, 
which modifies the PP in the HTS near the interface 
with a normal metal. Our results are consistent with the 
small imaginary component observed by STM (as a low 
transparency contact) as reported by Sharoni et al.flil] 
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FIG. 1: (color in online version) Theoreti- 

cal calculation(T=0):normalized conductance vs 
normalized voltage. (a) d PP, (100) contact, 
Z=0(black),0.1(red),0.4(green),0.5(blue)(b)d + is PP, (100) 
contact, z=0.5,fj=0(black),0.25(red),0.5(green),0.75(blue) 
|A| is the maximum amplitude of the PP. 



TABLE I: Summary of fit parameters and X ray data . 



No. 


Ca 


Xray" 


Tc b [K] 


Tc d c \K] 


Z 


A/ 


A/ 




1 





(100) 


91 


89 


0.34 


20.1 


16.3 


0.9 


2 





(110) 


86 


78 


0.39 


13.0 


10.7 


0.9 


3 


0.05 


(110) 


88 


77 


0.40 


16.3 


12.0 


1.7 


4 


0.2 


(110) 


76 


64 


0.46 


15.9 


8.5 


1.7 


5 


0.1 


(001) 


91 


82 


0.49 


17.1 


7.2 


1.8 


6 


0.05 


(001) 


90 


85 


0.52 


25.1 





1.4 


7 





(110) 


83 


63 


0.53 


11.4 


7.2 


3.2 


8 


0.1 


(110) 


91 


86 


0.60 


19.4 





2.1 


9 


0.1 


(100) 


79 


71 


0.68 


16.6 





2.8 


10 


0.1 


(001) 


85 


71 


0.68 


18.6 


0.8 


2.7 


11 


0.05 


(110) 


88 


77 


0.70 


18.1 


1.3 


6.1 



"(100) orientation was used in the fit for all contacts 
'Transition onset 
c Transition downset 
d meV 



In this study we have used Y\- X C , a x Ba2CuzOi -$ sput- 
tered thin films, with x=0, 0.05, 0.1 and 0.2. The growth 
procedure of these films were described in ref . [20I l2l| . 
The critical temperature of the films was determined 
by measuring the resistance versus temperature. X-ray 
diffraction was used to determine the global film orienta- 
tion (table 1). Sharoni et a/.jwjjj using STM topographic 
scans have shown that (110) films grown under the same 
conditions as those we have used, expose (100) facets. 
SEM and AFM pictures of the Ca doped (001) films show 
a-axis grains on the surface of the film. Thus, in all three 
film orientations (100), (110) and (001) (100) facets are 
exposed at the film surface. Ozawa et aL[15| have studied 
the electronic properties of the surfaces of (110) oriented 
YBCO films. They have found that the degradation time 
of (110) facets is significantly smaller than that of (100) 
facets. We therefore expect that in our films the chances 
of obtaining a good metallic contact with a (100) facet 
are considerably higher in comparison to a (110) facet. 



Our contacts were formed using a mechanically cut Au 
tip mounted on a differential screw. All measurements 
were taken at a temperature of 4.2K. 

We have analyzed the measured conductance spectra 
by fitting them to theoretical curves calculated using the 
formulas developed by Tanaka 5] . The procedure requires 
selection of the PP symmetry. We have tried the follow- 
ing possibilities: d, d+is, d+id X y and d+s. Within the 
selected symmetry the fitting algorithm has the following 
adjustable parameters : the amplitudes of the PP com- 
ponents, the barrier strength Z, the orientation of the 
surface in contact with the normal metal ((100), (110) 
etc.) and the contact's degree of directionality, namely 
the width of the tunnelling cone. The latter is dependent 
upon the barrier strength, Z. Following Wei et a/.^3|, in 
a high Z, tunnelling contact, we would expect the width 
of the cone to be around 20° while for a clean N/S con- 
tact we expect a value of almost 90°. Therefore in cal- 
culating the theoretical fitting curves, for our clean N/S 
contacts, we have used a 90° cone. The temperature for 
all calculated curves was set to 4.2K in accordance with 
the experimentally measured value. A life-time broad- 
ening parameter (Dvnesflij]') T was used to account for 
any smearing beyond the thermal one. We have found 
all our results to be best fitted using a d + is symme- 
try and (100) orientation, adjusting only the values of 
Ad, A s , Z and T(see table 1). Some of the curves can be 
fitted also using a d + id xy symmetry PP, however this 
requires using a narrow tunnelling cone in the fit, which 
is unreasonable for a metallic, low barrier, contact. As 
for the d+s PP symmetry, it requires using A s > Ad, 
thus, we find it less probable, as it fails to explain ex- 
perimental data obtained for high Z contacts jMllTj and 
would suggest that the s-wave channel is stronger than 
the d-wave one, which is in conflict with the findings of 
most experimental data [l8j |. 

In Fig. 2 we show three examples of our measured con- 
ductance data and the best fit curve for each of them. 
They are ordered by increasing contact's transparency, 
as determined by the Z value obtained from the best 
fit. The absence of a ZBCP suggests that for a d PP 
the data could be fitted only assuming a (100) oriented 
contact. Fig. 2a shows a maxima in the normalized con- 
ductance with a value of around 1.3 which one would 
expect could be fitted using the d symmetry, (see Fig. 
la). Indeed, the best fit parameters are given by Z=0.68, 
A = 18.6meFcos(20) + i0.8meV,T = 2.7meV and (100) 
orientation. In this case the additional is component 
needed to fit the data is very small and considering the 
smearing factor T = 2.7meV could very well be even zero. 
Fig. 2b shows two distinct features (marked by arrows 
in the figure), first a maxima at around ±7.4 mV, and 
second a distinct change of slope at around ±18 mV. It 
is impossible to reproduce this behavior for the case of 
a pure d PP, as it has only a single energy scale, but 
we can correlate these experimental features to the val- 
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FIG. 2: Normalized conductance vs. Voltage; T=4.2 K (circles) Theoretical fit (line). (a) Sample 10, Rn = 4057, A = 
18.6me^ cos(26>) + iO.SmeV, Z = 0.68, T = 2.7meV and (100) orientation(b)Sample 5,R N = 240, A = 17 .lmeV cos(26) + 
i7.2meV, Z = 0.49, F = 1.8meV and (100) orientation. Arrows point out to the manifestation of the two energy scales, maxima 
and change of slope(c)Sample 1,R N = 4.80, A = 20.1meV cos(26>) + H6.3meV, Z = 0.34, F = 0.9meV and (100) orientation . 



ues of A s and |A| = \J (A<j) 2 + (A s ) 2 respectively in the 
case of a d + is PP (see Fig. lb). For this contact, 
we find Z — 0.49, and the amplitude of the is compo- 
nent is enhanced. The other fit parameters are given by 
A = l7.1meVcos(26)+i7.2meV,T = l.SmeV, and (100) 
orientation, giving ^ w 0.4. Fig. 2c shows the con- 
ductance curve measured on our highest transparency 
contact, Z — 0.34. The amplitude of the maxima in 
the normalized conductance is around 1.7. Comparing 
to Fig. la, it obviously can not be fitted using the pure d 
PP (a peak amplitude of 1.7 can be reached using Z w 3 
but then the zero bias conductance would be lower than 
the high bias one). Using the d+ is PP we can repro- 
duce the higher peak amplitude without lowering the zero 
bias value (see Fig. lb). We obtain a best fit using 
A = 20.1mel/cos(26>) + il6.3meV, T = 0.9meV, and 
(100) orientation. As can be judged from Fig. lb the 
main effects of adding an is component are to increase 
the split between the conductance peaks and to enhance 
the maximum conductance. This trend is apparent in 
Fig. 2a trough Fig. 2c. 

In Fig. 3 we have plotted -j^j and j^j as a function of 
Z. As the samples have different critical temperatures we 
can expect variations in the values of A^ and A s even 
for a constant value of Z. To eliminate these changes we 
have normalized the two PP components, Ad and A s , by 
dividing them by |A|. The relative amplitude of the is 
component increases from values as low as 0-10% for Z 
around 0.7 up to values of around 60% for Z around 0.3. 
At the same time the relative amplitude of the d-wave 
component decreases. There is no obvious correlation 
with the level of Ca doping or with T c . The slope of ygr as 
a function of Z is maximum at Z 0.5. At low Z, the two 
normalized PP components appear to converge towards 
a value close to i.e. A s w A^. We thus see a change 
from a PP which is an almost pure d-wave for the high Z 
contacts to a PP with almost equal amplitudes of the two 
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FIG. 3: The relative weight of the PP components, 
■^■(circles) and m (squares) as a function of the barrier 
strength, Z. The low Z value to which the two normalized 
PP components appear to converge, (dotted line) 



components in the low Z regime, the crossover between 
the two regimes occurring around Z = 0.5. In Fig. 4 we 
have plotted |A| as a function of Tcdownset- The data 
is consistent with 2|A| = i]KbT c , where r\ = 6.0 ± 0.4. 
The fact that there is some scattering in the data is to be 




FIG. 4: Maximum amplitude of PP, |A| = ^/(A d ) 2 + (A s 
as a function of Tc downset (triangles) 
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expected as the films are never completely homogeneous. 
One must remember that |A| represents a local property 
of the film (at the contact), whereas T c measures a global 
property. 

We have shown that our data corresponds to a d+is 
symmetry PP, with the value of the imaginary compo- 
nent being a decreasing function of Z. A development of 
a subdominant imaginary PP (SIPP) near the surface of 
a d-wave SC was predicted theoretically by Tanuma et 
aL0 who preformed a self consistent calculation of the 
spatial dependence of the PP. They did not predict any 
change in the PP symmetry for a (100) contact. For 
the case of a (110) contact in which the predicted SIPP 
is maximal, the lower is the contact's transparency, the 
stronger is the reduction of the d-wave amplitude and the 
higher is the SIPP amplitude, i.e. the SIPP amplitude 
is predicted to be a decreasing function of the contact's 
transparency. This is inconsistent with our experimental 
findings. However, Tanuma et al. did not take into ac- 
count the possibility that an induced pairing amplitude 
can appear in the N side, i.e. a proximity effect, and the 
related effect this may have on the PP in the SC in the 
vicinity of the interface. The proximity effect between 
a HTS and a normal metal was studied by Y. Ohashi^ 
who predicted that a s-wave symmetry pairing amplitude 
is induced in N by the d SC. According to his prediction 
this should lead to a reduction of the d amplitude towards 
the N/S interface on the superconducting side. This ef- 
fect is maximized for a contact with the (100) face of a d 
SC and is enhanced in high transparency contacts. 

We believe that our findings can be explained as a re- 
sult of a proximity effect by making the following change 
in the proposal of Y.Ohashi. Indeed when a normal metal 
is in contact with a (100) boundary of a d-wave SC there 
appears a s-wave symmetry pairing amplitude in the nor- 
mal metal and a decrease of the d-wave amplitude in the 
HTS near the boundary. But in order to reduce the loss 
of condensation energy on the S side, we suggest that a 
is SIPP develops in the d-wave SC, in the vicinity of the 
barrier, in a way similar to that predicted by Tanuma et 
aL|7j for a (110) oriented contact. Considering our exper- 
imental findings, we conclude that YCaBaCuO has a sub- 
dominant s-wave pairing channel, otherwise the appear- 
ance of the is PP would have been energetically quite un- 
favorable. This subdominant channel exists for a broad 
range of Ca doping. Comparing our results to those of 
Gonnelli et al. [13], obtained on LaSrCuO, we find that 
though the symmetry used by Gonnelli et al. is a different 
one {d+id xy compared to the d+is) both experiments in- 
dicate that an imaginary component is added to the d PP. 
All of the spectra reported by Gonnelli et al. are for high 
transparency contacts (Z < 0.5) therefore no Z depen- 
dence of the two components can be inferred from their 
results. However the ratio of — T , w=rk w ~ 0.6 

^(A dxy ) 2 + (A d;c2 _ y2 )2 

that they find for doping values around the optimum level 



is not too far from the one we report here for the same 
doping levels in YCaBaCuO. It is important to point out 
that though we conclude that the large is component 
that we have measured for highly transparent contacts is 
a result of an unusual proximity effect, we do not suggest 
that this is the only scenario leading to the appearance 
of a SIPP in HTS. Various groups have measured a split 
in the ZBCP, explained by an appearance of a SIPP, in 
both planar [1J4 |22| and STM [1J| tunnelling junctions 
with (110) orientation, where the low transparency sug- 
gests that the mechanism leading to the appearance of 
the SIPP is unrelated to a proximity effect. A SIPP was 
also found by Farber and Deutscher|23j from measure- 
ments of the temperature dependence of the penetration 
depth in Ca over-doped YBCO samples, where no normal 
metal is involved. 
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